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Abstract 

We evaluate the Cold Nuclear Matter effects on J /iff production in pPb and PbPb collisions at the current LHC 
energy, taking into account the gluon shadowing and the nuclear absorption. We use the complete kinematics in the 
underlying 2 — > 2 partonic process, namely g + g — > J/if/ + g as expected from LO pQCD. The resulting shadowing is 
responsible for a large J ftp suppression in pPb and PbPb, and shows a strong rapidity dependence. 
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1. Introduction 

Relativistic nucleus-nucleus (AB) collisions are expected to produce a deconfined state of QCD matter - the Quark 
Gluon Plasma (QGP) - at high enough densities or temperatures. The J/ifr meson should be [1] sensitive to Hot and 
Dense Matter (HDM) effects, through processes like the colour Debye screening of the cc pair. A significant sup- 
pression of the J/ifj yield was observed at SPS energy by the NA50 experiment [2], and at RHIC by the PHENIX 
experiment in CuCu [3] and AuAu [4] collisions at y/s^^ = 200 GeV. The data recently taken at LHC in PbPb 
collisions at y/SMN = 2.76 TeV will provide results at a new energy scale, providing means to further test the available 
models. However, concurrent mechanisms - the Cold Nuclear Matter (CNM) effects - are known to already impact 
the J/ifi production in proton (deuteron)-nucleus (pA or dA) collisions, where the deconfinement can not be reached. 
Hence, the interpretation of the results obtained in AB collisions relies on a good understanding and a proper subtrac- 
tion of the CNM effects. Two CNM effects are of particular importance [5]: (i) the shadowing of the initial parton 
distributions (PDFs) due to the nuclear environment, and (ii) the breakup of cc pairs after multiple scatterings with the 
remnants of the incident nuclei, referred to as the nuclear absorption. In our previous works [6-9], we developed an 
exhaustive study of these effects. We confronted our results to the measurements from PHENIX [10] in <iAu collisions 
at -\/snn = 200 GeV, before giving our CNM effects estimates in CuCu and AuAu collisions. It is our purpose here 
to extend our results to pPb and PbPb collisions at the current LHC energy -\/saw = 2.76 TeV. 

As we have shown in earlier studies [6-9], considering the adequate Jfi/i partonic production mechanism - either 
via a 2 — > 1 or a 2 — > 2 process - affects both the way to compute the nuclear shadowing and its expected impact on 
the J /if/ production. From now on, we will refer to the former scenario as the intrinsic scheme, and to the latter as the 
extrinsic scheme. Most studies on the J/ift production in hadronic collisions are carried out in the intrinsic scheme. 
They rely on the assumption that the cc pair is produced by the fusion of two gluons carrying some intrinsic transverse 
momentum k T . The partonic process being a 2 — » 1 scattering, the sum of the gluon intrinsic k T is transferred to 
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the cc pair, thus to the J /iff since the soft hadronisation process does not alter significantly the kinematics. This is 
supported by the picture of the Colour Evaporation Model (CEM) at LO (see [11] and references therein) or of the 
Colour-Octet (CO) mechanism at a 2 s [12]. Thus, in such approaches, the transverse momentum Pj of the J /if/ entirely 
comes from the intrinsic k T of the initial gluons. However, the average value of k T is not expected to go much beyond 
~ 1 GeV. So this process is not sufficient to describe the Pj spectrum of quarkonia in hadron collisions [11]. 

In addition, recent theoretical works incorporating QCD corrections or s-channel cut contributions have empha- 
sized [13-15] that the Colour-Singlet (CS) mediated contributions are sufficient to describe the experimental data for 
hadroproduction of both charmonium and bottomonium systems without the need of CO contributions. Furthermore, 
recent works [16] focusing on production at <? + e~ colliders have posed stringent constraints on the size of CO con- 
tributions, which are the precise ones supporting a 2 — > 1 hadroproduction mechanism [11]. As a consequence, J /if/ 
production at low and mid Pj likely proceeds via a 2 — > 2 process, such as g + g — > J/if/ + g, instead of a 2 — » 1 
process 1 . This amounts to the bulk of the J /iff production cross section. Consequently, one is entitled to consider that 
the former 2 — » 2 kinematics i.e. the extrinsic scheme is the most appropriate to derive CNM effects at RHIC, and 
to provide predictions at LHC energy. In this work, we shall focus on the CNM effects expected at the current LHC 
energy Y^/wv = 2.76 TeV in the extrinsic scheme. The article is organized as follows: in section 2, we will describe 
our model and in section 3, we will present and discuss our results. 

2. Our approach 

To describe the J/if/ production in nuclear collisions, our Monte Carlo framework [6, 18] is based on the proba- 
bilistic Glauber model. The nucleon-nucleon inelastic cross section at xTsim - 2.76 TeV is taken to be ct^n - 
64 mb [19] and the maximum nucleon density to be po = 0.17 nucleons/fm . We also need to implement the partonic 
process for the cc production model that allows to describe the pp data, and the CNM effects. 

2.1. Partonic process for the cc production 

For Pj>2 - 3 GeV, most of the transverse momentum of the quarkonia should have an extrinsic origin, i.e. the 
J/iff's Pj would be balanced by the emission of a recoiling particle - a hard gluon - in the final state. The J /if/ would 
then be produced by gluon fusion in a 2 — > 2 process. This emission, which is anyhow mandatory to conserve C- 
parity, has a definite influence on the kinematics of the J /if/ production. Indeed, for a given J/ifi momentum (thus for 
fixed rapidity y and Pj), the processes discussed above, i.e. the intrinsic g + g — » cc — > Jfif/(+X) and the extrinsic 
g + g — > J/ifi + g, will proceed on the average from initial gluons with different Bjorken-jc. Therefore, they will be 
affected by different shadowing corrections. 

In the intrinsic scheme, the measurement of the J /if/ momentum in pp collisions completely fixes the longitudinal 

momentum fraction of the initial partons: x\^ = -^j= exp (+y) = x® 2 (y, Pj) with nij = ^M 1 + P 2 , M being the J /if/ 
mass. On the contrary, in the extrinsic scheme, the knowledge of the y and P T spectra is not sufficient to determine x\ 
and X2- Actually, the presence of a final-state gluon introduces further degrees of freedom, allowing several (x\ , X2) for 
a given set (y, Pj). The four-momentum conservation results in a complex expression of X2 as a function of {x\ , y, Pj)'. 



Equivalently, a similar expression can be written for X\ as a function of (x2,y, Pt)- Even if the kinematics determines 
the physical phase space, models are anyhow mandatory to compute the proper weighting of each kinematically 
allowed (x\,X2). This weight is simply the differential cross section at the partonic level times the gluon PDFs, 
i.e. g(xi,pF)g(x2,fiF)do-gg^j/,/, + g/dydPTdxidx2. In the present implementation of our code, we are able to use the 
partonic differential cross section computed from any theoretical approach. In this work, we shall use the Colour- 
Singlet Model (CSM) at LO at LHC energy, shown to be compatible [14, 20] with the magnitude of the P T -integrated 
cross-section as given by the PHENIX pp data [21], the CDF pp data [22] and the LHC pp data at y/s NN = 7 TeV. 



1 One may also go further and consider more than two particles in the final state, as expected from the real-emission contributions at NLO and 
NNLO [13]. It is clear from the yield polarisation [17] that these contributions start to dominate for Pt above 1 - 2m c . The effect of more partons 
in the final state is to increase the difference between the results obtained in both schemes. However the implementation of NLO and NNLO codes 
in a Glauber model with an inhomogeneous shadowing is not yet available. 
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2.2. Shadowing and nuclear absorption 

To obtain the J/iJ/ yield in pA and AA collisions, a shadowing-correction factor has to be applied to the J /iff yield 
obtained from the simple superposition of the equivalent number of pp collisions. This shadowing factor can be 
expressed in terms of the ratios Rf of the nuclear Parton Distribution Functions (nPDF) in a nucleon belonging to a 

nucleus A to the PDF in the free nucleon: Rf (x, Q ) = A ^ tldm(x Ql) , i = q,q,g . The numerical parameterisation 

of Rf(x, Q 2 ) is given for all parton flavours. Here, we restrict our study to gluons since, at high energy, the J /iff is 
essentially produced through gluon fusion [11]. Several shadowing parametrisations are available [23-26]. In the 
following, we shall restrict ourselves to EKS98 [24], which is very close to the mean in the current evaluation of 
the uncertainty [26] on the gluon nPDF and exhibits a moderate antishadowing. We postpone the propagation of the 
uncertainty on the gluon nPDF to the CNM effects evaluated at LHC energy for future studies. 

The second CNM effect that we are going to take into account concerns the nuclear absorption. In the framework 
of the probabilistic Glauber model, this effect is usually parametrised by introducing an effective absorption cross 
section <x abs . It reflects the break-up of correlated cc pairs due to inelastic scattering with the remaining nucleons from 
the incident cold nuclei. The value of cr abs is unknown at LHC. At high energy, the heavy state in the projectile should 
undergo a coherent scattering off the nucleons of the target nucleus [27], in contrast with the incoherent, longitudinally 
ordered scattering that takes place at low energies. As argued in [28, 29], this should lead to a decrease of <x abs with 
increasing Y*aw- The systematic study of many experimental data indicate that cr abs appears either constant [30] or 
decreasing [31] with energy. Hence, we can consider our estimates [8, 9] of cr abs at RHIC energy as upper bounds for 
the value of cr abs at LHC. We choose three values of <x abs that should span that interval (<x abs = 0, 1.5,2.8 mb). 



3. Results and discussion 



In the following, we present our results for the J /if/ nuclear modification factor due to CNM effects in the extrinsic 



= 2.76 TeV: R AB = dN J '% /(N con )dN J J ', where dN J '*(dN J ' ¥ ) is the 
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sheme in pPb and PbPb collisions at ~\[t 
observed J/ifr yield in AB = pPb, PbPb (pp) collisions and {N co n) is the average number of nucleon-nucleon collisions 
occurring in one pPb or PbPb collision. Without nuclear effects, Rab should equal unity. 
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Figure 1: (Color online) J/ifr nuclear modification factor versus y in pPb and PbPb collisions at s/s^n = 2.76 TeV, 
using EKS98 [24] gluon shadowing parametrisation and three values of cr aft , (from top to bottom: 0, 1.5,2.8 mb) in 
the extrinsic scheme. For PbPb collisions, the y-dependence is shown for various centrality selections. 

In Fig. la, we show R pP \, versus y. The curve with no absorption allows to highlight the strong rapidity dependence 
of the shadowing. We can also notice that the shadowing alone should already be responsible for a quite large amount 
of J I \f/ suppression, up to 34% at y — 4. This is expected due to the very small x -region in the gluon nPDF that 
becomes accessible at LHC energy (down to 10~ 5 ). Fig. lb shows that the y-dependence of 7?p b p b is similar for all the 
centrality bins, with a dip at mid-y. This shape is the opposite of the one obtained at RHIC energy [6, 8], with a peak 
at mid-y. Here, ^p b p b is systematically smaller at mid-y than at forward-y. This is also illustrated on Fig. 2, with the 
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Figure 2: (Color online) J/if/ nuclear modification factor, /?pbpb, in PbPb collisions at y/s^N — 2.76 TeV versus Af palt , 
using EKS98 [24] gluon shadowing parametrisation and three values of the nuclear absorption cross section in the 
extrinsic scheme. Rp\,pb is shown for two different experimental acceptances in rapidity, \y\ < 0.9 and \y\ e [2.5,4]. 



centrality dependence of RpbPb for two regions in y. This behaviour of the CNM effects may partially - or completely 
- compensate the opposite effect expected from cc recombination, with a maximum enhancement at y — 0. Overall, 
one may observe a RpbPb rather independent of y resulting of two y-dependent effects. 
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